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ABSTRACT 

An orbit perturbation procedure is applied to the description of 
monochromatic, larKC-amplltude, electrostatic plasma wave propagation. 

In the lowest order approximation, untrapped electrons are assumed to 
follow constant-velocity orbits and trapped electrons are assumed to 
execute simple harmonic motion. The deviations of these orbits from the 
actual orbits are regarded as perturbations. The nonlinear damping rate 
and frequency shift are then obtained in terras of simple functions. The 
results are in good agreement with previous less approximate analyses. 
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1. INTRODUCTION 


There is currently strong interest in the nonlinenr behavior of 
plasmas Involving large-amplitude waves, or high levels of turbulence. 

The relevant theory has advanced from analyses of wave-wave and wave- 
particle interactions bused on second- or tliird-order perturbation aeries 
in the wave amplitude (Kadomtsev 196'^; Sagdeev and Galeev I 969 ; Tsytovich 
1970 ), to the inclusion of higher-order nonlinear effects, which become 
important as the level of perturbation increases. The theory is quickly 
limited by mathematical difficulties, one of which stems from time secu- 
laricies in the perturbation series caused by trapped particles (Davidson 
1972 ). Since the larger the wave amplitude, the more particles are trapped, 
the series rapidly becomes non-con vergent . Furthermore, trapped particles 
oscillate with periods pi’oportional to a fractional power of the wave 
amplitude (proportional to the square root of the amplitude for particles 
at the bottom of a monochromatic wave potential). Thus, when trapped 
particle dynamics are significant, the series does not give a proper 
approximation. 

A number of theories have been developed to include trapping effects 


for a coherent, electrostatic wave, when the wave amplitude is not too 

1/2 

large, i.e. when a[ --- (v /vJ(ep)/T) ] « 1 lor a one-dimensional 

cp t 

Maxwellian electron velocity distribution; here v and v [= (T/ra) ] 

Cp ^ 

are the wave phase velocity and electron thermal velocity, respectively 

(O’Neil I 96 D; Al'tshul and Karpman I 966 ; Tanluti I 969 1 Imamura, Suglhara 

and Taniuti 1969; Dee and Schmidt 1970; Bailey and Denavit I 97 O; Oei 

and Swanson 1972; Suglhara and Kamimura I 972 ). The various theories 

have predicted different time evolutions of the wave depending on the 

value of Y , where Y is the linear Landau damping rate and 

L ^ B L 


n r = (ekE /m)^‘^‘^l is the electron bounce frequency. When 7 /n 3 ; 
3^* ' 0 L 13 

the wave damps out according to the Landau (l9^W:j) treatment, before 


trapped electrons can execute any significant bounce motions. For larger 
wave amplitudes, the trapped electrons tend to flatten the velocity 


distribution locally, so that the damping rate decreases with time from 

its initial value, Y . For a wave of even larger amplitude (typically 

y /oj„ « 1 ), the initial Landau damping rate decreases rapidly at about 
Li B 
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a- n , and osclllatcg thoreaftor as tho trapped electrons bounc'j. This 
oscillation eventually disappears^ and the wave evolves towards a steady- 
state Qernstein-Greene-Kruskal mode (19!>T)* wave frequency in this 

state is different from tho linear value, tho difference being proportional 
to Djj (Mnnheimer and Flynn 1971 » Morales and O'Neil 1972; Dewar 1972 I 
Lee and Pocobelli 197? and I 97 J) , Tho foregoing characteristics are in 
general agreement with experiments (Malmberg and Wharton I 96 (' ; Franklin, 
Hambergor and Smith 1972; Vidmar, Malmberg and Starke 197 !j); account 
being taken of the fact that the theories treat temporal rather than 
spatial variations, and with computer simulations (Armstrong 19671 Tsai 
197 ^! Canosa and Gazdag 197^M Matsuda and Crawford 1979)* 

There is considerable interest in cases lor wliich the condition 
o: « 1 is unnecessarily restrictive. For example, sideband instabilities 
are prominent ior 0 ; > 1 (Wliarton, Malmberg and O'Neil I 968 ; Franklin, 
Hamberger, Ikozi, Lampis and Smith 1972; Jahns and Van Hoven 1973) • In 
fact, a « 1 is more restrictive than tho small-amplitude condition of 

p 1/2 

the usual perturbation theory for the non-resonant region, (e E“/nT) « 1 

O /p U 0 

or a « (v (Bud'ko, Karpraan and Shklyar 1972). 

For a 1 , the wave may have quite different characteristics from 

those described for a « 1 . For example, consider an arbitrarily large 

wave amplitude, so that Y « 1 even .for low phase velocity. We might 

anticipate that this wave will propagate at an arbitrarily low phase 

velocity, and reach a steady state by trapping. However, this does not 

happen; since the distribution function increases rapidly with decreasing 

velocity, the number of trapped electron increases rapidly with the wave 

amplitude. These electrons take more energy from the wove, and thus 

enhance the wove damping (Armstrong 19^7: Dawson and Shanny I 968 ; Sato, 

Ikezi, Takahashi and Yamashita I 969 ; Nakamura and Ito I 97 I 1 Tsai 197^ i 

Vidmar, Malmberg and Starke I 97 ? i Sugihara and Yajnanaka 1975 » Canosa 1975)* 

Consequently, a wave of amplitude larger than a certain value damps out, 

the precise value depending on the phase velocity. Dawson and Shanny 

( 1968 ) have shown that this eventual total damping occurs for all values 

of amplitude when v /v^ ^ 3 . 5 * 

CP t 


2 


It iB extremely difficult to extend riKorous theories for q « 1 
to the case of a > 1 . Here, we shall sacrifice exactness, and apply 
a perturbation procedure to the case of a « 1 . Specifically, we 
consider a coherent, monochromatic, electrostatic wave, and separate the 
dynamics of untrapped and trapped electrons, foliowiriK O'Neil In 

the lowest approximation, we assume the untrappod and trapped electron 
orbits to bo of constant velocity and simple hannonic motion, respectively. 
We then consider the difference between these oi'bits and the actual orbits 
as perturbations. 

In § 2 , the damping rate and frequency shift are obtained from the 
lowest-order orbit approximation. In §3; on appropriate perturbation 
theory is developed. This is applied in §’l and §') to untrapped and 
trapped electrons, respectively, to include the effect of the orbit 
perturbations. In §6, our results are compared with those of Morales 
and O'Neil (1972)> ond some general conclusions are drawn. 
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LOWEST OUDEH APPROXIMATION FOR NONLINEAR DISPERSION RELATION 

Consider an infinite^ hoinoKcnooiis, collisionless plasma, with a 
neutralizing background of immobile positive ions, driven by a single, 
longitudinal wave E(x, t) ^ shall restrict our interest 

to an initial value problem, wliero k is a real constant, and will assume 
that temporal variations of the wave amplitude, E^ , and frequency (hence 
the wave phase velocity) are small. In the wave frame, the electron 
velocity distribution function in phase space, f(x, v,t) , is described 
by the '^lasov equation 


»i„ „x MiL^LSi , „ . (1) 

3t 5x III 0 ov ' ' 


Together with (l), Poisson's equation 

O 00 X/P 


'0 


_ p 


s/-/ 


^Af(x, V, t) oxp(ikx) , (2) 


constitutes the complete set of equations. Here Af(x,v, t) = f(x, v, t) - fQ(v), 
f^(v) is the initial distribution function, and X( - 2it/k) is the wave- 
length. Equation (1) implies conservation of the distribution function 
along the electron orbit. 


X V , V ^ — E sin kx . (3) 

m u 

Although f is expressed conveniently in ( Lagrangian) coordinates 
moving with the electrons according to (3), is better described by 

fixed (Eulerian) coordinates. This difficulty may be minimized by finding 
an invariant of the electron motion, and using this as an independent 
coordinate. If E^ is infinitely small, we may neglect the right-hand 
side of the v-equation in (3)^ regard v as a constant of the motion, 
and expand (1) in terms of E^ to the desired order. This is the pro- 
cedure of the usual small-signal theory (Sogdeev and Galeev 1969)* For 
large amplitude, we can find several invariants, e.g. energy, the trapping 
parameter (O’Neil 1965)> ihe action variable (Goldstein I95O), and 
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ti-nnsform the (x, v) -coordlun toH to this new coordinnto for (1) nnd (?) . 
ThlB had boon done by nuiny outliors cltod In §1. Tho ti’nnsfonnntlon Is 
rather Involved and, in most cases, tho inteurntlons in (?) can not be 
performed analytically. In order to simplify tho procedure, we first 
obtain the orbit equation from ( j) !, y multiplyini; tho v equation by v , 
and integratinK with respect to time, 



H sin^ 



(>0 


where H is the total energy of an electron, anti v^_( = Uj^/k) is the 
trapping velocity. The electron orbit described by (It) is drawn in 
figure 1(a) . 

To lowest-ordor, we may approximate ()|) by a constant velocity 
orbit for tho untrnppod electrons (ll > mv^) . We can then follow the 
procedure of the usual small-signal theory to solve the Vlasov equation. 
Neglecting higher-order mode-mode coupling, we obtain the well known, 
linearized solution of (?), 






[cos k(x-vt) - cos kx] 




where the prime denotes the velocity derivative. For trapped electrons 

p 

(H < mv^) , we replace the wave potential by a parabola as shown in 
figure 1(b). Near the bottom of tlie potential well. It is quite accurate. 

The largest difference arises near the wave crest. However, electrons 
travel very slowly near the crest and do not exert a significant effect 
on the wave. Therefore, we consider this potential difference a perturbation. 
In the zeroth order approximation, (*i) reduces to 

+ (kx)^ = 

which describes a circular orbit in the (kx, v/v^) -plane. 
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with those npprovimntloiiH, tho Vlnsov criuntioii niny bo written ns 


ai 






2 { sin kx T \ 

■5’‘ 1-iS- - 57 ■ 


(7) 


The chaructorlstic curve of the lolt-hnnd side is the clrculnr orbit of 
( 6 ). The right-hand side represents tho per tiirbntlon duo to tho difference 
between tho simisoldnl and parabolic potentials. Changing tho variables 
(x^v) to clrculnr coordinates (r, iji) defined by 

lex r sin ij; , v/v^_ r cos v(i , ( 8 ) 


( 7 ) can bo put in the fonii 

bp 3 f / sln(r sltnj;) 

"51 ‘ n ^ ‘ b y r 



r cos'jf 



/ 


(9) 


whore F{r,i;(,t) = f(x( r,-j/, t), v(r, (jf, t), t) . In the zeroth order approxi- 
mation for trapped electrons, we neglect tho right-hand side of ( 9 ). Tho 
solution is tlien 

F(r,.j,,t) - • (10) 


Following Morales and O'Neil (I972), wo assume that the wave ampli- 
tude is small (a « 1 ) so that, in the trapped region, ^q(v) can be 
expanded about v = 0 os 


f^(v) = f^(0) + vf'(O) + ^ f"(0) 


( 11 ) 


Substituting this into (lo) gives 

AF(r,t,t) s F(r,^,t) - F(r,^,0) 


= 2 rv 


.f'(0) 


sin 


(¥) “■(' - ¥) 




fj,^( 0 ) sin a)jjt sin( 2 tir-a)jjt) 

( 12 ) 
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In ol'dor to use (' ) iitui (L ), wo linvo to HopiU’iito cont rilni ttons to 
(l') of im trapped mul Lnippod oiocLroiis, K^, ■ . I’or tivi iin trapped 

electron contribution, we dotennine tlie lower Imiind | vj - v,^ of the 
urtrnijpcd rcKion l^y considoi’lng tlie ni'oa between the x-nxls and the 
orbit of the electron durinj; one wave period. This ni'eu is tlie action 
variable (Gold.stcin 1') o). I'rom (^t), wo see tluit the minimum of tlii.s 
urea occurs when il mv‘^_ and takes the value in tlie (kx, v/v^,)-plnne. 
Since the untrnppod electron orbit is of constant velocity in the lowest 
order approximation, the lower bound is j^iven by k>v^j/v^ B , or 
Vj^j (l^/J^)v . The untrapped electron contribution to (p) is then 



-u 


Substituting (‘j) gives as 

The untrnppod electron contribution, , tlioroforo becomes 



B 



where the barred integral sign denotes the Couchy principal value. In 
the real part, the first integration of the cosine term phase-mixes to 
zerc. For the second integration, wo may use (11), In the imaginary 
part, most of the contribut:,on comes from near |v| « v^ , and we may again 
use (ll) for 1 q(v) • Carrying out the integration reduces (I5) to 




sin kVyt V 


- ^"^0 (^f ("• ’ I ^ 


(16) 


where Si (a) is defined by 


si(a) = 


a 


(17) 


Changing variables (x, v) to (r, iji) gives the trapped electron 
contribution to (&) as 


nvu 

= -2iek 


‘hJ 


rdr Fj^(r, t) , 


( 18 ) 


n 

j(r,t) = J AF(r,ijf,t) exp(ir slnijf) , 


2 

where the upper bound of the r-integratlon is determined by na /2 = 8 
1/2 

or a = h/K . Substituting (12), and performing the ijr -integration 
(Cradshteyn and Ryzhik, I 965 ) yields 




r f \ 

- 


r^J. 


r^(r)(l-coB 2oijjt) + lv^f^(0)rj^(r)sin . 
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Llsln(' the (lerivntivo l’oniiul»i for Bessel I'unctions (Abramowitz and Stegun 


the r-lntogrations in (I 8 ) can bo performed to obtain 

E„ . 2^ j,(n) (^f ^ V 

' ^ ' (ai) 

The Poisson equation involves the sum of (K^) and (21). Substitution 
yields the dispersion relation 


e = e + Oe = 0 , 

L * 


(22) 


where e is the linear dispersion relation; e and 6 e ore given by 
L b 






6e 


i‘M •fi 

-(r) 




, / sin kv t\ 3 

i (1 - "k"yt ■) " TT 


-t- i It 






1 - - Si(kv,,t) + J^(a) sin u t 

It ' U ' It 2 B 


. ( 23 ) 


We now let d = ;Uq + 61 . -t- iV , where ji^ satisfies the linear 
dispersion relation (®^ " ^ ''xpand (22) about . To first 

order in 6 s ^ 6 ;u ^ and Y , we obtain 

5 / sin(1.27.>yt)\ 

0 -^ = -^5 (1 - . 0.5, 


;) (1 - COS 2«)gt) , 


^ = 1 - 0.6)+ Si ( 1,27 'Ojjt) + 1.32 sin aigt , 
L 


(21+) 
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p // 

wliorc n 5 (:ii /k)*^v f (o) (&e , the Laiuiau dampinp: coefficient 

0 pro 1, 0 -j /p 

Y = - jt(.a /krfA(0)(ae /biO"^ , and ) = 0.41 and J (4/a = O.lT 

Ij P 0 L U c; J 

linve been substituted. The first bracketed terms for 6;j> , and the first 

two terms for Y , ore the untrapped electron contributions; the remaining 

terms are the trapped electron contributions. 

The frequency shift and damping coefficient of (24) are shown in 

figure 2. They exhibit the essential features of experimental observations 

(Malmberg and Wharton I 967 ; Vidmnr, Malmberg and Starke 19T5) : 

Y/Y decrease within or near i> t « ; oscillatory behavior follows with 
Li n 

periods of about 2,u„t for 6 i,/fj and n t for Y/Y . The time- 

u 0 13 Lt 

asymptotic value, ® reasonable agreement with the 

more accurate value of 4 d/q^ = - 1.63 obtained by Morales and O'Neil 

(19T2). 

To ^ain insight into the relative importance of trapped and untrapped 
electron contributions, they are shown separately in figure 2, as well as 
combined. The sinusoidal behavior of the former is well understood from 
energy and momentum considerations for the trapped electrons and the wave 
(Morales and O'Neil 1972). The behavior ox the latter has not been 
clarified previously: In the linear theory, the constant velocity orbit 

is assumed for tlie entire region of velocity. In the present case, of a 
lai’ge-araplitude wave, it can be assumed only for | v| > . The Poisson 

equation therefore assumes the form of a finite Pourier transform [see (I 3 )]. 
This contributes the (l - sin x/x) and Si(x) [x = 1.27 terms in 

(24). The initial increase of 4 o/fi^ and Y/Y^^ of the trapped electron 
part is thus offset by the contribution from the untrapped electrons. 

It should be noted that, whereas Morales and O'Neil show ^u/q^ =: 0 
and Y/Y - 1 for some time before they decrease, figure 2 shows 
^ ‘j/Qq > 0 and Y/Y^ > 1 for n^t ^ 0,6n . This difference is mainly 
due to our approximation that all trapped electrons have the same bounce 
frequency, and hence exchange energy and momentum with the wave coherently. 

The actual bounce frequency is, of course, lower for electrons closer to 
the separatrix. If we include this different bounce frequency, the first 
peaks of 6j,/f^ and Y/Y will be reduced; the amplitudes of their 
oscillations will be decreased; the maxima and minima will be delayed. 
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r/x, 






and the oscillations will eventually decay away, ns demonstrated by 
Morales and O’Neil (197*0. To include those effects, we shall take 
account of the difference between the actual and the lowest order orbit, 
as a perturbation effect, in succeeding sections. 
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I 


I 
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3. ORBIT PERTURBATION 

So far, we have considered tlio lowest order orbits: constant velo- 

city for untrapped electrons, and cii’cular orbits for trapped electrons. 
Under the influence of the wave, however, an untrapped electron changes 
its velocity as it travels. Similarly, a trapped electron changes its 
radius and the angular frequency because of the difference between the 
parabolic and actual sinusoidal wave potential. Wo shall assume the 
time-scales of these changes to be large compared to the time for an 
untrapped electron to traverse one wavelength, or for a trapped electron 
to complete one cycle of its trapped orbit. 

We write the Vlasov equation in oper.ttor form, 

9f(x, t) a 

+ Lf(x,t) = 0 , L = ^ (25) 

whore x represents both spatial and velocity coordinates, g is the 
generalized velocity, and the summation convention on the repeated index 
is understood. We divide f into a slowly-varying part, (f^ > Jind a 
rapidly-varying part, f^ , wliere ( ) denotes a time average over the 
shorter of the two time-scales. We assume (g) = 0 . [If this is not 
the case, we need only consider x to be a coordinate moving with (g) , 
and (g) then disappears from (25)]. Taking the time average of (25) 
gives the secular part. Subtracting this from (25) gives the rapidly- 
varying part, 

df 

<Lfj_) = 0 , + (l-A)Lf^ + L<f) = 0 , (26) 

where 1 and A are identity and averaging operators. These equations 
are identical to those used in the renormalization procedure of plasma 
turbulence (Dupree I966 ; Orszag and Kraichman 19^7 : Frisch I968; 
Weinstock I 969 ; Rudakov and Tsytovich 1971; Tsytovich 1972; Benford and 
Thomson 1972." Misguich and Balcscu 1975)* 

Following Tsytovich (I 972 ), we introduce an average collision 
operator, v , such that 


15 


(27) 


v<f) s - <LI^) , 


and rewrite ( 26 ) as 


a<f> 


= v<f> . 


(28) 


This expression is then integrated to obtain 


(f> = <f(0)) + / dt^ v(t^)<f(t^ 

■{) 


)> 


t t ^ 

= /i +/* dti v(t^^) +y* dt^ y 

' *r\ r\ 


dtg v(t^) v(tg) + ... 


0 


4- t 1 

t 1 n-1 


/ d^ 2 *** f v( ^]^) v( , , , v( tj^) + ''*1 (l(0)) . 

V\ ‘'n ’ 


0 0 


(29) 


Unless indicated otherwise, throughout this section the argument of 
f is time. Introducing a time-ordered exponential function defined by 
(Kubo 1962 ) 

t ^ t t 

exp^y dt^v(t^) = ^ ^y* dt^...y dt^0[v(ti)***'^(^n)l 

0 n=0 0 0 

^n-1 


y dt^... / dt^ ^ti)...v(t^) , 


( 30 ) 


we may cast ( 29 ) in the compact form 

<f(t)> = U(t,0) <f(0)> , 


( 31 ) 
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whore the average propagator, U , is given by 

(32) 

the operator O denotes the chronological ordering in v^s t ®nd U 
satisfies the group property U(t,tQ) = u(t, t^) U(t^, t^^) (Weinstock lSK>9) • 
For the rapidly-varying part, f^ , we may use the qiiasilinear form 
of (26) (Benford and Thomson 1972: Misguich and Balescu 197!?) 

Sf 

vfj_ L<f) = 0 . (33) 

The solution is then 

t 

fj^(t) = - j U(t,t^) L(t^) U(t^,0) <f(0)> , (34) 

■q 

where we have assumed tlie initial condition f^(0) = 0 • Substituting 
this into ( 27 ) gives 

t 

v(f) = / (i'(^) ^(^>^j^) ^(^^^^) • (35) 

0 

The effect of the modified orbit on f^ is contained in U , which 
can be found by iterating (32) and (31?)* U = I, (34) reduces to the 

result for the lowest order orbit. 


U( t, t^) 


t 

= "«Po / 


dt. 


At,) 


“0 
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UNTKAPPKl) EU-XTKON CONTRIBUTION TO MJCCTUIC I'lEU) 

In 0 coorcltnotc system moving with the lowest ortlor straight lino 
orbit (y^w) , clofined by 

y x~vt , w : V , ( 36 ) 

the Vlasov equation (l) takes the fonn of (2'') with 

- I Kq k(yl-wt) (t - y . (37) 

To carry out tho iteration of (32) and (3‘0, we first let U = I . Noting 
that sin k(y+wt) and ( tS/3y - 5/t)w) commute, wo con write (31^) “s 


t 

'' ('^ 57 - s )/ "-1 

0 


t^ ^sin k(y+wt) sin k(y+wt^)) |r “ 






i) a 


3y 


\ I dT cos kw 
ow 


(('-•') 57 - 17) 


( 38 ) 


Since the time-scale of slow variations is assumed to be much longer 
than 1/kw , the integrand contributes only for a very short time. Tho 
upper limit of the time integration may therefore be taken as t -* <n , 
Noting that 


/ 


dt cos kwt = jtii(kw) 


. 0 , 


CO 



sin kwt 


0 


kw ^ 


(39) 


which follow from the fact that kw r 0 for untrapped electrons, we can 
reduce ( 38 ) to 


V 



(kw) 



1 a 


)• 


(40) 
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Next; we may obtoin U(t, t^) from (32) aa 



where the ttnio-ordei‘ing subscript, 0 , has been dropped because the 
exponent no longer involves the tiine-intograt j on. Wo now substitute (41) 
in (34) to obtain . For the U(t,0) operating on (f(0)) , wo first 
note that, since (f(0)) is homogeneous, 5/by ^ 0 , Therefore, 

Ll(t^,0) -= I . For the same reason, !•( ^j^) becomes -(e E^/m) sin k(y-twt^) ( 5/5w) . 
In the first U , we hove 5/5y = ± ikw . Transforming back to the (x, v)- 
coordinate consequently reduces (34) to 



Since we assume that the time for on un trapped electron to traverse 
one wavelength, T 1/kv , is much shorter than other variations, we 
may put T » T in (42) • Therefore, we neglect the second term in the 
square bracket in (42). The first exponential function in the integrand 
then takes the form exp (A-i-b) , which con be expanded as (Misguich and 
Balescu 1974; Louisell 19^4) 

exp(AfS) = exp(A) exp(B) exp(- [A, B]/2) ... , (43) 

where [A, B] = AB - BA . Since [ A, B] = ®(®q) > hence negligible, we 

may factorize the exponential function. Since kvT » 1 , the most 
significant contribution to b/bv comes from operating on exp(ikvT) . 
Therefore, we replace b/5v by ikT . The Fourier transfonn of (42) 
then becomes 


19 


II -III 



(410 

Since the I’irat teitn in the exponent gives only n slow vnrlotlon, the 
inteqrutlon yioids upproximotoly 



Oo) 


p 

Note that, in the limit j)„t/(v/v ) . o(lf) - 0 , we recover iront 
(4|j) the linear result of (14). For simplicity, we shall neglect this 
small term in the denominator. Since we see that, when (4|)) is substituted 
in ( 13 ), the dominant contribution of tlie exponential term is from the 
small-v region, we shall write v^ for v in the first term of the 
exponent in (4^). With those considerations, (lb) is finally modified to 


CO 



( 46 ) 
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TItAPPED ELECTRON CONTUIUUTION TO ELECTRIC KIEU) 


In order to cost tne Vlnaov oquntioti for trnppocl oloctrons, ( 9 ), in 
the fonn of ( 2 ^)), wo nvornKo the rinlit-hniui side of ( 9 ) with respect to, 
, over the ron«o (-«, .i), nnd subtract tliis nvorn( 40 d port from both 
sides. We then obtain 


&F 5f 

5t ^ 


+ 


hF 0 , 



r 




(>‘7) 


whore the functio»is R and S ore defined by. 


n . 

R(r, 'jf) -s ^ |sin 2 ^ cos i)( sin (r sin iit)| , 


S( 


( 


cos S'!! I- — sin i}t sin (r sin '{() 




(48) 


and is Riven by 


Bl 







( 49 ) 


Trapped electrons traveling according to 

; = , r - 0 : ^ = :i^ 3 _t -I- -j,y , 

feel the rapidly-varying force exerted by tlie L^term in (47). To fiti' 
the orbit correction due to this force, we shall assume in 

( 90 ) for simplicity. In terms of the coordinates > ^^^7) can 

then be reduced to the foiir, of (2‘ ), with L given by 


r ; r 


0 ' 


(:^o) 


/I V 





(:^1) 


where R aed S are given by (1^8) with (r,\it) replaced by 
according to ( 50 )* 
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Substltutinj; ( 1} inul U I in (i ), nnd LnltiiiK the time nvcriiKC 
over ono bounco poriotl, jj > obtain 


V j(lT ^L(t) L(t-T)> - j, V ^ , 

I 


r..? 


6 ^ 

2 

1 a 

!ir &'r 

r 

1,' 


+ — 
)• 

J (r ) ~ ^ J, (2r ) 
(r t ' r 1 J 



0 0 


t' 

(J 


Trunsfoiming 


to 


('^7) 

now hocorncjs 




5t ' 

Sf 

- (t- 

1- u)f , 

(--3) 


where Q - w (i -i- 6- -t- 6 ) niul 0^ is ^iven by ('.2) with r 
replaced by r . 

The dominant orbit correction if= obtninod by nOKiectiiiK the right- 
hand side of (V3), yioldinf; the equntion 


3f 

'ST '■ 6 ;;. 


0 , 


whoso solution is 






rv^_ cos (^if 







Substituting this in (l8) gives 


Fko(r,t) iv^£'(v) rj^(r) sin n^^t - vV(o) 

(%) 

For the quasi linear correction, we write ('3) as 


dF^ 

^ ^13 "Sr ■ ■ ''‘''n • 


(‘■7) 
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AssuminB &Qj^/br 0 ^ for simplicity, we obtain 


t 

0 

X cos(<)(H 2 jjT)j , ( 58 ) 

Substitutins this in (I 8 ), and usiiiB (11), (48), and integration formulas, 
(Crarishtcyn and Hyzhik gives 


F„.(r,t) = - i 


IJ 

Vriy(O) Jd6 jrJ^(r)cos(n,jt-29) + 2 J^( r) t-P( r, 9) cos(njjt-Q) j 


■«">(%:)/' 


d9 r^j^(r) sin29 - sii.(20jjt-9) 


2rjj^(r)j2(r) siii2njjt H- 2i’Q(r,9)cos(nj^t-9)j . 


Here, p(r, 0) and Q( r, 9) are given by 


P(r,9) - JQ^2r Bin |j - ^2r cos 


Q(r,9) = ^2r sin |) cos(n^t - (2r cos |) pin (n^t - 1 1 


Using standard Bessel function identities (Abrnmowitz and Stegun I 965 ) we 
may express P and Q in series form as 


L 


23 


1 


P(r,9) . COS (2nH-l)9 , 


n=0 


“ 1 

Q(r,G) = 2^^ J2(n-)-l)^"’^ sin[npt (2nn-l)Q]- J2^(r)8in[oj3t-(2in-l)0j 


n =0 


(61) 


The terms with ii ^ 0 can be shown to make neglifjlbly small contributlonSi 
Retaining n=0 terms in (A. 5 ), substituting in (A. 3 ), and performing 
0 -integrntlons gives 


F^(r,b) , lv^,i;(0) Bln n^t 


2 


wliere, since P, , is small, we liave used ;n„ « n„ • Adding this to 

^ kl ' B B 

Fj^O (‘-^*^) gives 

^ ^ko-" "^kl 


iv/'( 0 ) rj^(r)[l-i-e^(r)j sin n^t 


- vV(0) 


(63) 

The contribution to the electric field due to trapped electrons, 

E , is then obtained from (I 8 ) as 


= - 2 iE„ I 7 ^ 


co 2 /* 

=0 (ir) ■■^Ji(--)[l+\(r)]Bin n^t 

0 


(61^) 


Eq /'Ji. \2 
~ \k 


2 /* 


2h 


1 


1 
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6. CONCLUSIONS 

The results of §l^ and §5 are the equations (46) and ( 58 ) • 
them gives the dispersion relation (22) with given by 


Combining 


* 

fte = f^f 


sln(l.2T rujjt) 
1.2T iDjjt 


) 


- 0.5 

0 

\2 


a 

>Jdr r^|jg(r)[l+e^(r)] - J^( r)6^( r)j (1-cos 2n^t) 


ia(^) f'(0)|exp[^- (0.39 V^) 


O.ShJ dr r^ J^(r) [l+6^(r)]sin n^tj 


W 


The frequency shift and damping rate are found by expanding e about 


The frequency shift 

•0 = COq . 

We obtain 

flg) _ 

- 2.55 

^0 


/ r 2l to t) \ 

= - 2.55 - exp[- (0.39 a>3t)2J j 

a 

O.^J dr 


^ = exp[^- ( 0.39 aigt)^] - 0.64 Si (1.27 to^t)! 


0.64/ dr r^j^(r)[l+6^(r)]sin CgC . 

0 (65) 
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Figure 3 shows 6oi/n„ oncl Y/Y calculated from { 6 ^^), with contri- 

0 L 

butious of trapped and untropped electrons shown separately and combined. 

Comparison with the results of Morales and O’Neil (1972) generally shows 

good agreement. The asymptotic value of Coj/Qq = “ 1.62 is very 

close to their result of = - 1.63 • figure indicates that 

the untrapped electron contributions to and y/Y change rapidly 

for o)jjt remain substantially constant for ai^t > a . As was 

mentioned in §2, the origin of this characteristic can be easily .^eon by 

dividing the velocity integration over the region |''[ < v^ , as shown 

in (I 9 ) the first integral, jver the entire velocity region, gives the 

linear dispersion relation and the Landau damping rate; it is the second 

integral, over | v| < > that gives the time variation of the untrapped 

parts of aud Y/Y . In the work of Morales and O'Neil, the 

0 L 

second integral has been combined with that of trapped electrons in (18), 
and is consequently observed in the subtraction of the linear charge 
density from that of the actual perturbed density. 

The damping rate, Y/Y , and frequency shift, , are proportional 

1/ ^ ^ 

to ^0^^^ ■’ • Moi'^les and O'Neil explain this from 

momentum and energy considerations. If four conjugate points keep their 
symmetry during motion in phase space about the potential well (which is 
the case of a wave with constant amplitude and phase velocity), the total 
momentum change experienced by them derives from the odd-order velocity 
derivative terms in the Taylor series of (11) ; the contributions from the 
even order derivative terms cancel out because the momentum is odd in 
velocity. On the other hand, since the energy is even in velocity, the 
total kinetic energy change of these conjugate points derives from the 
even-order derivative terms. Since we keep only terms up to the second 
velocity derivative of the distribution function, Y/Y and 6 oj/q are 
consequently proportional to ^q(^) ^ larger amplitude 

wave than has been considered here, it is apparent from the above picture 
that, because of additional contributions from higher-order terms of the 
Taylor series, the frequency shift and damping rate can be larger than 
the present results. 
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FIG. 3 (a). Frequency shift obtained from the orbit 
perturbation. Dots are results of Morales and 
O'Neil (1972). 
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Wave amplitude and phase velocity chanKes may also enhance the 
frequency shift and damping rate. It should be remembered that the 
present theory assumes constant wave amplitude and phase velocity, i.e, 
the momentum and energy gain by those four conjugate points during 
acceleration is returned to the wave during deceleration. However, when 
the wave damps initially, some of the trapped electrons which have gained 
momentum and energy during initial acceleration are detrapped, and a 
corresponding amount of momentum and energy is not recovered by the wave 
during the next decelerating period. Similarly, the phase-velocity 
change resulting from the frequency shift disrupts the symmertrical 
motion of the four conjugate points. The damping rate and frequency 
shift are then no longer necessarily functions of only odd and even 
orders, respectively, of the velocity derivatives of the distribution 
function. 

Although it has been used for a monochromatic wave, the perturbation 
procedure developed in §3 is identical to that previously used lor 
random waves. This is an important point since the literature developed 
in §3 contains many analyses specialized to one or other of these two 
extremes, but there is no easy way to combine them. The method used 
here effectively links them, and is potentially applicable to such 
important situations as slight randomness accompanying coherence, and 
vice versa. 
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